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Abstract 

If the cosmological dark matter has a component made of small primordial black holes, they 
may have a significant impact on the physics of the first stars and on the subsequent formation 
of massive black holes. Primordial black holes would be adiabatically contracted into these 
stars and then would sink to the stellar center by dynamical friction, creating a larger black 
hole which may quickly swallow the whole star. If these primordial black holes are heavier than 
~ 10^^ g, the first stars would likely live only for a very short time and would not contribute 
much to the reionization of the universe. They would instead become 10 - 10^ A/q black holes 
which (depending on subsequent accretion) could serve as seeds for the super-massive black 
holes seen at high redshifts as well as those inside galaxies today. 
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1 Introduction 



The first stars in the Universe mark the end of the cosmic dark ages, reionize the Universe, and 
provide the enriched gas required for later stellar generations. They may also be important as 
precursors to black holes (BHs) that coalesce and power bright early quasars. The first stars are 
thought to form inside Dark Matter (DM) halos of mass 10^ Mq at redshifts z = 10 - 50 

(Abel et al. (2002); Bromm et al. (2002); Yoshida et al. (2003)). These halos consist of 85% DM 
and 15% baryons in the form of metal-free gas made of H and He. Theoretical calculations indicate 
that the baryonic matter cools and collapses via H2 cooling (Peebles & Dicke (1968); Matsuda et al. 
(1971); HoUenbach & McKee (1979); Tegmark et al. (1997)) into a single small protostar (Omukai 
& Nishi (1998)) at the center of the halo (for reviews see Ripamonti & Abel (2005), Barkana & Loeb 
(2001) and Bromm & Larson (2004)). 

It is interesting to study the effects on the evolution of the first stars due to the large reservoir 
of DM within which these stars form. The first protostars and stars are particularly good sites for 
this investigation because they form inside the highest density environment (compared to today's 
stars): they form at high redshifts (density scales as (1 + z)^) and in the high density centers of 
DM haloes. Previously, two of us (Spolyar et al.(2008)) studied the effects of Weakly Interacting 
Massive Particles on the first stars: we found that the annihilation of these particles could provide 
a heat source for the star that stops the collapse of the protostar as well as potentially dominates 
over any fusion luminosity for a long time. 

In this paper, we consider instead the effects on these first stars of a different candidate for the 
DM: Primordial Black Holes (PBHs). These are small black holes that may be formed in the very 
early universe (see the next section for more detail) and may exist in sufficient abundance to provide 
the DM seen in the Universe today. The masses of PBHs that explain the entirety of the DM range 
from (lO^"^ - 10^6) g; heavier PBHs up to 1 Mq may provide still an interesting fraction of the DM. 

We discuss the implications that PBH DM would have on the physics of the first stars, the so 
called Population HI stars. These stars could range from ^ (1 - few 100) Mq. First, we compare 
various possible heat sources due to PBHs with the ordinary heat from stellar fusion of the stars. For 
the properties of the Pop III stars, we use results computed by Heger & Woosley. Specifically, for a 
100 Mq (10 Mq) Pop. HI star, we take the central temperature to be 1.2 x 10^ K (9.6 x 10^ K), the 
central density 31 g/cia^ (226 g/cm^), the radius 7 Rq (1.2 Rq), and the stellar fusion luminosity 
to be 



We find that the ordinary stellar fusion luminosity dominates over the heat sources due to PBHs, 

which include accretion onto the BHs, Hawking radiation, and the Schwinger mechanism. 

Instead, we find the interesting result that PBHs inside the first stars may sink to the center and 
form a single BH, which may accrete very rapidly and swallow the whole star. The phenomenon 
is relevant for PBHs heavier than about 10^^ g, because the corresponding timescale for dynami- 
cal friction turns out to be shorter than the typical stellar lifetime, while it is less interesting or 
completely negligible for lighter BHs. So, for Mpbh ^ 10^^ g, the lifetimes of Pop. Ill stars may 
be shortened, with implications for rcionization of the Universe as well as for the first supernovac. 
In addition, since the stars are inside much larger haloes, they can in principle accrete even more 
matter (depending on the accretion mechanism). Thus, the end-products of the scenario are BHs of 
masses 10 — 10''' Af©. These may be the seeds which produced the super massive BHs seen at high 
redshifts; the Intermediate Mass Black Holes; as well as the black holes at the center of every normal 
galaxy today and whose origin is as yet uncertain. Possible mechanisms of production of superheavy 
BHs are reviewed in Dokuchacv et al. (2007). In addition, although the PBH swallowing the star 
shortens the star's lifetime and its contribution to reionization, the newly formed hole can become 
a new, alternative source of ionizing photons. 

The rest of the paper is organized as follows. In Section II, we briefly review the physics of 
PBHs, that is, how they can be formed in the early Universe and what current constraints on their 
cosmological abundance are. In Section III, we discuss the behavior of individual PBHs: how many 



= 6.5 X lO^^erg/s (100 Mq), 
= 4.2 X lO^^erg/s (10 Mq). 
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of them are expected inside a single star (via adiabatic contraction), what is the luminosity due to 
accretion onto the PBHs, and what is the timescale for their size to double. We also investigate 
alternative mechanisms for generating luminosity by these small PBHs. Then in Section IV and V 
we turn to the most important part of the paper. We study the dynamical friction that pulls all the 
BHs into a single larger BH at the center of the star, and then watch this single large BH accrete 
the entire star surrounding it on a fairly rapid timescale. We conclude with a discussion in Section 
VI. Throughout the paper, we use units with c = 1. 

2 Physics of Primordial Black Holes 
2.1 Production mechanisms 

PBHs may be formed in the early universe by many processes (Zcldovich & Novikov (1966), Hawking 
(1971), Carr & Hawking (1974), Crawford & Schramm (1982), Hawking (1989), Polnarev & Zem- 
bowicz (1991), Dolgov & Silk (1993), Jedamzik (1997), Rubin et al. (2000), Dolgov et al. (2008)). 
For a general review, see e.g. Carr (2003). The earliest mechanism for BH production can be fluc- 
tuations in the space-time metric at the Planck epoch. Large number of PBHs can also be produced 
by nonlinear density fluctuations due to inhomogeneous baryogenesis at small scales (Dolgov & Silk 
(1993), Dolgov et al. (2008)). If within some region of space density fluctuations arc large, so that 
the gravitational force overcomes the pressure, we can expect the whole region to collapse and form 
a BH. In the early Universe, generically, BHs of the horizon size are formed, although it is also 
possible to form much smaller BHs (Polnarev & Zcmbowicz (1991), Hawking (1989)). BHs can also 
be produced in first and second order phase transitions in the early Universe (Crawford & Schramm 
(1982), Jedamzik (1997)). Gravitational collapse of cosmic string loops (Polnarev & Zembowicz 
(1991), Hawking (1989)) and closed domain walls (Rubin et al. (2000)) can also yield BHs. The 
masses of PBHs formed in the above mentioned processes range roughly from Mpi (BHs formed at 
the Planck epoch) to Mq (BHs formed at the QCD phase transition). 

The basic picture is that energy perturbations of order one stopped expanding and recoUapscd 
as soon as they crossed into the horizon (Zeldovich & Novikov (1966), Hawking (1971), Carr & 
Hawking (1974)). The maximal mass of PBHs is set by the total mass within the cosmological 
horizon, i.e. Mhor = M^i/A"^ at any given energy scale A at which the BH forms. This is also the 
expected mass scale of a BH in most early Universe scenarios for the production of PBHs (it can be 
at most a factor of 10~^ smaller (Hawke & Stewart (2002))). Thus 



where we assumed a radiation dominated Universe, with the effective number of relativistic 
degrees of freedom and Tf the temperature of the Universe at time tf. 

2.2 Observational constraints 

PBHs in the mass range Mpbh ~ 10^^ — 10^^ g can be good DM candidates. A number of constraints 
restrict the mass to this range. PBHs with an initial mass smaller than about 5 • 10^^ g are expected 
to be already evaporated due to Hawking radiation; moreover their presence in the early Universe 
can be constrained by observations for Mbh ^ 10^ g (lifetime r > 1 s) (Novikov et al. (1979)). For 
MpBH ~ 10^^ g, there are strong bounds as well, at the level of flpsH ^ 10~^, from the observed 
intensity of the diffuse gamma ray background (Page & Hawking (1976)), so they may be at most 
a tiny fraction of the non-relativistic matter in the Universe. For larger masses, constraints can be 
deduced from microlensing techniques (Alcock et al. (2000); Tisserand et al. (2007)) and dynamical 
arguments (Carr & Sakellariadou (1999)), which exclude the possibility that the whole cosmological 
DM is made of BHs heavier than 10^^ g, even if they still may be an important component. For 
example, the PBH to DM mass ratio in the Galactic Halo would be smaller than 0.04 for PBHs in 
the mass range 10^° - 10^^ g and than 0.1 for the mass range 10^'^ - 10^^ g (Tisserand et al. (2007)). 




(3) 
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On the other hand, for the mass range 10^^ — 10^^ g, there are currently no clear observational 
methods of detection. For Mpbh ^ lO" - 10^" g, the presence of PBHs can be inferred from the 
femptolensing of gamma ray hurts (Gould (1992), Nemiroff & Gould (1995), Marani et al. (1999)), 
but the constraint is weak, roughly Qpbh ^ 0.2; in addition it holds only for uniformly distributed 
DM and is not easy to extend to the more realistic case of clumped DM. The same mass range might 
be covered by future gravitational wave space antennas, from the gravitational interaction of PBHs 
with test masses of the laser interferometer (Seto & Cooray (2004)), but the expected detection 
rate for LISA is too low and only a further generation of space detectors might put non-trivial 
constraints. According to recent work Abramowicz et al. (2008), the PBH mass range 10^^ — 10^^ g 
remains unexplored and thus allowed. However, further constraints raise the lower bound to roughly 
IQie _ 1017 g (Bambi et al. (2008a)). 

We present results for PBHs with mass Mbh = 10^"* g but show the scaling for other PBH 
masses in the 10^^ — 10^^ range. Our results are qualitatively the same for PBHs of any mass in 
the allowed range. For heavier PBHs up to e.g. IMq, the results will be somewhat different and 
discussed in the discussion section. 



3 Primordial Black Holes inside the Star 

In this section we study the behavior of the PBHs inside the star. We estimate the total mass in 
these objects, as well as the luminosity and timescale for accretion onto individual PBHs. 

3.1 Total Mass in PBHs inside the star 

The first stars form at the centers of IO^Mq DM haloes. As a starting point we assume an initial 
Navarro-Frenk-White profile (Navarro et al. (1997)) for both DM (85% of the mass) and baryons 
(15% of the mass). As the gas collapses to form a star, it gravitationally pulls the DM (in this case 
PBHs) with it. We use adiabatic contraction (Sellwood & McGaugh (2005)) to find the resultant 
dark matter profile inside the star (Spolyar et al. (2008)) 

PDM ~ 5 {nt cm3)"-8 GeV/cm^ , (4) 

which is independent of the nature of DlvQ. Here, ni, is the mean baryon density inside the star. 
It should be noted that adiabatic contraction is not a relaxation process. Instead as the baryons 
collapse to form a star, they gravitationally pull the DM will them, so that the DM density inside 
the star increases. Hence, the DM evolves on the timescale of the baryons. Ideally, instead of using 
the adiabatic approximation, it would be desirable to run an N-body simulation. At present this 
is technically not possible. Regardless, adiabatic contraction should give a reasonable approxima- 
tion and is widely used formalisrrH. In addition, the formal requirements to apply the adiabatic 
approximation hold during most of the evolution of the baryons. For a mean baryon number density 
nt, Ri lO^'' cm^'^, the DM to baryon matter mass ratio of a typical Pop. III.l star is at the level of 
lO^"'. The number of PBHs inside the star is roughly 



, Mbh J VlOO Mq, 

where is the mass of the Pop. HI star. More precisely (modeling the star as an n = 3 polytrope), 
we find for a 100 Mq (10 Mq) star that the total mass in PBHs is 

Mtot,PBH — 6.3 X 10^" g (100 Mq) , (6) 



^This is the result of a calculation for DM density in the first stars that we performed with WIMP dark matter in 
mind, but exactly the same result holds for any type of DM including BHs which are orders of magnitude larger. 

^Our original work on adiabatic contraction in the first stars was performed using a very simple assumption of 
circular orbits only. However, in follow-up work, two of us participated in a paper (Freese et al. (2009) ) in which 
we performed an exact calculation including radial orbits. The results changed by less than a factor of two, so that 
we feel comfortable using Eq. (|4]l. In that same paper we also considered a core alternative to an NFW profile as 
our starting point for the adiabatic contraction and, again, obtained essentially the same result. Our results for DM 
densities in the first stars appear to be quite robust. 
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Mtot,PBH = 

4.1xl0^«g (10 Mq). (7) 
3.2 Accretion onto the PBHs from stellar material 

In this paper we study the effects of PBHs on the stars on the main sequence, once they have fusion 
proceeding in their cores. The PBH effects are much more important during this stage than during 
the protostellar collapse phase. Since they are surrounded by a high density stellar environment, 
they accrete and emit radiation. As a maximum possible value, the accretion luminosity for a single 
PBH cannot exceed the Eddington limit 

Le = — = 6.5 . 10 j erg/s , (8) 

where aru is the Thomson cross section. is the luminosity at which the outwards radiation 
pressure compensates the gravitational attraction and stops the accretion process. Clearly, the 
Eddington luminosity is proportionate to mass. In this case, the mass has been conservatively set to 
the mass of the BH (Mbh)- In fact, the mass should include the optically thick gas surrounding the 
BH. Under this restriction, the maximum stellar luminosity from PBHs inside one star is realized 
when the accretion luminosity of every BH is at the Eddington limit, i.e. 



M, 

iooIv/q 



'^E,tot = NbhLe - 10^^ ( ) erg/s. (9) 



Since LE.tot oc MpBH,tot, the upper bound on the power emitted by PBHs is independent of the 
BH mass. This accretion powered luminosity is at least a few orders of magnitude smaller than the 
expected stellar luminosity for Pop. Ill stars, 4- lO'^^ erg/s (6- lO'^^ erg/s) for 10 and 100 Mq (Freese 
et al. (2008)) stars respectively. The extra heat produced by accretion onto the PBHs inside the 
star has thus a negligible impact on the physics of the star. 

As the PBHs accrete more matter and become more massive, the Eddington limit increases and 
the BH accretion luminosity becomes more and more important in the energy balance of the star. 
The Bondi accretion rate is (Bondi (1952)) 

MbhV fl^e^V'^ f Pb 



M.=4.i.|p...l.4.10-^^j ^^jg/s. (10) 

Here Rb = 2GnMbh /v'^ is the Bondi radius. The quantity v is the typical velocity of the particles 
of the accreting gas with respect to the BH, and should account for both the thermal velocity of 
the particles Vp = ^'ST/nip, where T is the local temperature of the star, as well as the BH orbital 
velocity vbh = \J GiqM^ir) /r, where M*(r) is the stellar mass within a distance r from the center. 
Close to the center v k, Vp, but for large r this relation is no longer true; instead, the BH orbital 
velocity may reduce the accretion rate, even by an order of magnitude. We here take v — Vp and 
use the Bondi formula to find an upper limit on the accretion rate, recognizing that this value may 
well overestimate the true accretion rateH. The differential equation Mbh = ck-^'-^b// has solution 

Mb nit) - , , , (11) 



1 - aMot 



where Mq is the BH mass at t = and 



aMo = 1.6 . 10-1^ ( ^) ( ( , ) s-i (12) 

\1024 gy \ T ) \l g/cm3 ' ^ ' 



^Moreover, the Bondi formula holds in the ideal case of perfect spherical symmetry. In realistic situations, the 
non— zero angular momentum of the accreting gas and the presence of other effects (magnetic fields, turbulences, etc.) 
may diminish the accretion rate, since La must be smaller than Le- The case of accretion onto BHs is however a 
complex phenomenon, because BHs have an event horizon and in principle may be capable of swallowing an arbitrary 
amount of matter without exceeding the Eddington luminosity (Begelman (1978), Begelman et al. (2008)). We will 
take the Bondi accretion as an upper limit (Begelman (1978)). 
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is the inverse of the characteristic accretion time of the BH. The accretion time scale is thus not 
shorter than 



It is possible for even a single PBH with Mbh > 10^'* g inside the star to eat the entire star. 
Such a case was discussed in Begelman (1978) in the context of a super-massive star capturing a 
BH in a bound orbit. The current scenario differs due to the fact that we are interested in the role 
of PBHs on Pop. Ill stars and their effects on cosmology (e.g. reionization); here the PBHs are 
thought to comprise at least some measurable fraction of the DM in the universe and are therefore 
present in the haloes containing the Pop. Ill stars before these even form. If the PBHs do not 
comprise the entire DM, then the PBH mass could be larger than we have discussed heretofore, 
though contributing only a small fraction of the critical density. 

As we will show below, the maximal accretion rate computed here is somewhat slower than the 
rate for the formation of a larger BH at the center of the star; all the effects combined thus lead to 
a big BH at the center. 



3.3 Other mechanisms for energy release by PBHs 

One may be also concerned about two other mechanisms in which PBHs can release energy: Hawking 
radiation (Hawking (1975,1976)) and positron emission (Bambi et al. (2008b)). 



3.3.1 Havirking radiation 

The luminosity due to Hawking radiation is maximal for the smallest mass BHs. We thus consider 
the (unrealistic) possibility that all the cosmological DM is made of PBHs with mass Mbh = 10^"* g. 
The Hawking luminosity per BH from 7, and fj,^ emission is 7- 10^* erg/s (Page (1976)) and their 
total contribution to the power of a 10 Mq star would be at the level of 10^^ erg/s, roughly 2 or 3 
orders of magnitude smaller than the ordinary stellar luminosity, 4- 10^^ erg/s. If the mass of the star 
were 100 Mq, the relative contribution would be smaller, because the stellar luminosity increases by 
a factor 100, while the BH luminosity increases by a factor 10. Higher Hawking luminosity would 
demand smaller PBHs. However, if the PBHs had an initial mass Mbh = 10^'^ g, their lifetime 
would be T < 10^ yr, that is much shorter than the age of the Universe when first stars formed. 
Thus fusion luminosity always dominates over the Hawking radiation. 



3.3.2 Schwinger effect 

The second mechanism, positron production due to Schwinger effect at the BH horizon, has been 
recently discussed in Bambi et al. (2008b). Because protons are much more massive than electrons, 
it is much easier for BH to capture protons. Whereas the protons fall right into the BH, the electrons 
interact frequently via Compton scattering on their way into the BH and are prevented from falling 
in as easily. Hence the BH builds up a positive electric charge. For a BH mass AIbh < 10^" g, the 
electric field at the BH horizon can exceed the critical value for vacuum stability, i.e. Ec = rn^/e, 
so that electron-positron pairs can be efficiently produced (Schwinger effect). Then, electrons are 
back-captured while positrons escape. The net result is to convert protons of the surrounding plasma 
into 150 MeV positrons. The accretion rate of protons is (Bambi et al. (2008b)) 

„,...(^)'(-)-(_»^).. 

We note that mechanism is not the same as Bondi accretion, and that the expression above is not 
obtained from eq. PU)) . By contrast, Bondi accretion is the accretion of gas where particles collide 
with one another other, losing their tangential velocity but gaining radial velocity towards the star. 
This hydrodynamic approximation is applicable if the characteristic length scale is larger than the 
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mean free path of particles. Here, the size of the BH is smaller than the proton mean free path, Ap, 
and we consider protons at distances r < Xp with small velocities, so they are gravitationally bound 
to the BH. These protons lose energy by bremsstrahlung or synhrotron radiation near the BH and in 
this sense they are not non-interacting. The picture is very much diffent from the hydrodynamical 
one and the calculations of the proton accretion rate can be found in Bambi et al. (2008b). Once 
equilibrium is reached between the accretion rate and the Schwinger discharge rate, the luminosity 
per BH is roughly (Bambi et al. (2008b)) 

For Mbh = 10^*^ g, this equation would then imply that the total Schwinger luminosity is roughly 
10^^ erg/s for a star of mass 10 - 100 Mq, which is comparable to the fusion luminosity for 10 Mq 
stars given in Eq. ([2]) but far below the fusion luminosity for 100 Mq stars given in Eq. ([T]). However 
this value of the Schwinger luminosity is never reached, because the rate for proton capture is 
10^^ s^i (6 X 10^° s"^) for a 100 Mq (10 Mq) star, while the rate to create the e+/e~ pairs is the 
product of the production rate per unit volume, ~ m^, and the volume of the region around the BH 
in which the electric field exceeds the critical value Ec- The latter is a spherical shell of thickness 
about 1/me, so the volume turns out to be r^/me, where = 2GnMbh is the BH gravitational 
radius. The pair production rate is ^ nilrg ~ 5 x 10^^ s~^ for a 10^° g BH (the Schwinger discharge is 
fastest for this BH mass) . Hence the equilibrium between the capture and the Schwinger mechanism 
is reached for a Schwinger luminosity that is several orders of magnitude lower than given above for 
the stellar mass = (10 — 100) Mq. Thus fusion always dominates over the Schwinger efi^ect as a 
heat source. 



4 Formation of a larger black hole at the center of the star 
via Dynamical Friction 

4.1 Main Sequence Star 

The most important phenomenon associated with the PBHs inside the first stars is the formation 
of a larger BH at the center. It is well known that gravitational interactions cause every heavy 
body moving into a gas of lighter particles to lose energy by dynamical friction (Binney & Tremaine 
(1987)). Thus, PBHs inside a star are expected to sink to the center of the star, eventually forming 
one single large BH. 

We will use Chandrasekhar's dynamical friction formula to compute the timescale for the PBHs 
to sink to the center of the star. If we assume that the gas of light particles has a Maxwellian velocity 
distribution with dispersion a, then the deceleration of a BH moving at a velocity vbh with respect 
to the rest frame of the fluid is 



dt 



Vbh = -47r Gj^Mbh Pb In A 



VBH 



^BH 



evi{X) 



2X exp(-A:2 



(16) 



where X = v b h / {'sf^cr) , erf is the error functioUjOt, is the density of the background particles and 
In A w hi{M^/MBH) is the Coulomb logarithm Q There are two possible regimes, depending on 
whether vbh is larger or smaller than the velocity dispersion inside the star, a = ^jTjm (Binney & 



'The actual definition of Coulomb logarithm is (see Binney &: Tremaine (1987)) 

h -2 

InA = In 



where bmax is the maximum impact parameter, a is the mean square velocity of the gas and m the molecular 
weight. Numerical simulations show that bmax can be assumed of order the orbital radius of the object, say _R. Since 
~ GffMt(R)/R, a reasonable estimate of A is Mi,{R)/Mgjj. 
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Trcmaine (1987)). Here T is the local gas temperature and m is the molecular weight. The factor 
in the square brackets 

n^..)^erf(X)^ ^^-P(-^^) (17) 

tends to unity for vbh 3> cr and tends to v^Qj^l2\p2/na^ for vbh ^ 

The vector equation (|16p can be rewritten as the following two scalar equations: 

r = 5 ^ ^ - -fivBH)r, (18) 

j = -i{vbh)J, (19) 

where r is the distance of BH form the star center, J ^ r^O \s the BH angular momentum per unit 
mass, 9 is the azimuth angle, vbh — + J"^ jr^ ^ 

M,{r) - / dhytir) (20) 

is the stellar mass inside radius r, and 

^9,, , . erf(X) - 2Xexp(-A:2)/y^ ^ ^ 



"BH 



Since the characteristic gravitational time scale 



1/2 / 3 \ 

1900 (ilZS^) s (22) 



is much shorter than the lower limit on the characteristic dynamical friction time scale 




"--^ - ^.GlMBBPblnK - ^ ■ ' ^ M ^T^T^. ) I ^7" 1 ( TTT ) ^ , (23) 

we can approximately solve eqs. ([TSl [T9|) as follow^. We may neglect the last term in the r.h.s. 
of eq. and assume approximate equality ~ GjvAf,(r)r. Using this result we can integrate 
eq. (|19p. which now takes the form: 

VBH = 5 VBH : (24) 

Vbh^^df 

and calculate the time of capture of small BHs at the stellar center. The result depends upon the 
initial velocity of the BH which we may estimate assuming that the BH is on a circular orbit of radius 
r determined by the stellar mass M*(r) interior to radius r, i.e., vbh = \/GN^I^r)/r. We find that, 
in the outer regions of the star, vbh ^ cr, in which case eq. ([24]) scales as vbh — —(t^/{tdfv^ij). 
In the inner regions of the star, near the stellar center, we find the opposite limit of vbh ^ c, in 
which case eq. ((24|) scales as vbh = —vbh / {'^V^tttdf) instead. Thus, in the latter case, the time 
of BH formation is about 

Tf « 




^The reader might also be concerned whether we can neglect the third term in eq. JTSj when considering the 
opposite limit as vbh goes to zero. In this case, a goes to vbh in eq. II23I I. Again, the third term in eq. I I18I I is 
completely negligible and even more so than when eq. II23I I depended upon a. 
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where is the initial PBH velocity, so « a and implies Rm ^ lO^*' cm, while u;^^ is the final 
PBH velocity, when Rf — 4 - 10^ cm, that is, when the orbit of the PBH is equal to the Schwarzschild 
radius of the final BH. In the case vbh ^ c, the timescale becomes 




1 



(26) 



which can be quite a bit longer than the one for the case vbh ^ f for v'j^j^ moderately larger than 
v'^^. As shown later, this is not a problem, because we always have a sufficient number of PBHs at 
small radii, where vbh ^ c. 

The case of very eccentric orbits does not significantly change the picture. A simple estimate can 
be obtained assuming radial motion and constant matter density pb- In the absence of dynamical 
friction, the motion of the PBHs can be treated as an harmonic oscillator with period Tg and velocity 
~ {Ro/tq) cos{t/Tg), where Ro is the maximum distance from the center of the star and t is the time. 
Since the maximum velocity exceeds 3 • 10^ cm/s for Rq > R.t ^ 10^^ cm, the equation of motion of 
PBHs inside the radius i?* is basically 

T T 

(27) 



TDF 

The differential equation is that of an underdamped harmonic oscillator: 

r{t) - e"*/^^"'^ cos{Lut + S) , (28) 

where w ~ ^/Tg. We find t/ — 2t£)fIii (^"^), a timescale which is actually shorter than in the 

circular case. Thus we expect that the result in eq. ([25|) is a reasonable estimate of the timescale. 

To obtain more accurate quantitative estimates of the dynamical friction timescale on which the 
PBHs sink to the center of the star, we did numerical calculations assuming that the star can be 
modeled as an n = 3 polytrope, which is known to roughly reproduce the stellar properties of a 
star dominated by radiation pressure. We can then obtain density and temperature profiles for a 
100 Mq star which are plotted in fig. ([1]). If one does the full stellar structure of a star of a Pop. HI 
star, the exact answer is different than found assuming a polytrope. The difference is on the order 
of at most tens of a percent. 

Subsequently, we can now compute the timescale for the case of Mbh = 10^^ g; luckily. The 
resultant timescale can easily be scaled to other BH masses since it is inversely proportional to 
Mbh- To be specific, we investigated the case of a 100 M© star. We found that the transition 
from fast to slow BH velocity (relative to gas particle velocity) takes place at i?c ~ 2 • 10^" cm. 
As explained above, the dynamical friction for BH outside of this radius is proportional to 1/vbhj 
while, for smaller radii it is proportional to vbh- Roughly 50% of the BHs are initially inside the 
radius Ri = 1.4 • 10^^ cm; these BH take 1 • 10** yr or less to sink to Re- (We have also computed 
the timescales for infall for BH coming in from different initial radii Ri to the same value of Re, our 
results are shown in Table ([T]). Subsequently each BH takes another ~ 5 • 10^ yr to sink from Rc to 
Rf = A ■ 10^ cm. The latter is the Schwarzschild radius of the final BH at the center of star. Thus 
the timescale for half of the PBHs to form a single large BH at the center of the 100 Mq star is 
roughly 

Thus for Mbh > 10^^ g, in a 100 Mq star, the timescale for the formation of a large central BH 
is less than a million years, which can have a significant impact on the evolution of the star. We 
note that, once the central BH mass is ~ lO^^g, the (fastest possible) accretion timescale in Eq. 
becomes comparable to the dynamical friction timescale Eq. 1^^; the result of both effects is a single 
large BH inside the star. 

If the mass of the star is 10 Mq, the sinking time is not significantly different. 
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Ri (cm) 


M,{Ri)/AU{R,) 


p{r,) (g/cm^) 


Time (yrs) 


4.7-10^^ 


1.0 


0.5 


205,000 


2.4-10" 


0.9 


3.1 


30,000 


1.4-10" 


0.5 


8.7 


11,000 


6.5 - lO^" 


0.1 


17.3 


5,000 



Table 1: Numerical results of the timescale for BHs to move from a variety of initial radii Ri to a 
smaller radius Rc ^ 2 ■ 10^" cm. ]\'h{Ri) / M^{R^) is the fraction of the mass of the star inside the 
radius which is equal to the initial fraction of PBHs inside the radius Ri. p{Ri) is the stellar 
density at Ri. The mass of the BH has been fixed to lO'^'^g 



Additional PBHs from outside the star may also fall onto the central BH via dynamical friction. 
For a baryon density profile that scales as Pbir) oc r~^ '^ outside the star, we find that the dynamical 
friction timescale is 

where is the initial radius of the infalling PBH and this equation has been computed in the fast 
BH regime with vbh oc v~g\j. Thus it takes a very long time for dynamical friction to be effective at 
pulling in BH from typical radii in the minihalo. From closer in, the timescale can be significantly 
shorter, e.g., it takes 150,000 years for a 10^'' g PBH to go from 3 x 10^^ cm (~ 10 times the radius of 
the star) to the center. However, the amount of mass in PBHs inside this radius is 2.4 x 10^^ g, more 
than an order of magnitude less than the amount already in the star from Eq. ([5]), and is therefore 
negligible. Thus dynamical friction does not pull in significant mass in PBHs from outside the star. 



4.2 Protostellar Phase 

One may ask whether dynamical friction is already effective during the protostellar phase, long before 
the Pop. HI star comes to exist on the main sequence. Early on, there is a collapsing molecular 
cloud which is very diffuse and becomes more and more dense as it cools via molecular hydrogen 
cooling. The protostellar clouds stop collapsing once they become protostar nuggets with 10~^ Mq 
in mass, hydrogen densities of 10^^ cm~^, and radii ~ 5 x 10^^ cm (Yoshida et al. (2008)). In the 
standard picture of Pop. HI star formation, there is then accretion onto these nuggets until the stars 
reach ^ 100 Mq and go onto the main sequence. 

Can the PBHs already sink to the center during this earlier phase and cause the protostar to go 
directly to a BH, avoiding the main sequence phase altogether? Inside the protostar, the appropriate 
regime for dynamical friction is that of slow BH, with vbh oc vbh- Such protostellar clouds have 
much lower densities than the subsequent Pop. Ill stars on the main sequence, and consequently 
are ineffective at causing the PBHs to slow down via dynamical friction. We have checked that the 
timescale is simply too long for PBHs to play any role during the collapse of the protostellar cloud, 
unless the PBHs are much more massive than have been considered in this paper. However, once 
the nugget forms, the baryon density is high enough to trap PBHs of mass > 10^^ g with dynamical 
friction. At this point the Kelvin- Helmholtz time ~ T£)i? ~ 10 yr. The amount of DM (PBHs) 
inside the nugget is ~ 10^^ g, so that the initial central BH is only this big. However, it quickly eats 
the rest of the 10"^ Mq protostar, and presumably can grow at least to the value of the original 
1000 Mq Jeans mass of unstable material. 



5 Eating Pop. Ill stars 

We have shown that the PBHs can sink to the center of the star and form a single larger BH in 
a reasonable timescale (for MpsH > 10^^ g) to change the evolution of the star. We now need to 
address the subsequent fate of the star: can the BH really accrete at the Bondi rate and swallow the 
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whole star quickly? Alternatively, does the radiation pressure from the accretion luminosity slow 
down the accretion rate and make the star have a normal evolution? In general, the accretion of 
matter onto an object with a solid surface (e.g. a neutron star) is limited by the radiation produced 
by the accreting gas. 

La = r?M (31) 

where rj is basically the gravitational potential per unit mass on the surface of the object. Neverthe- 
less, in the case of accretion onto BH, the picture is more complex and the phenomenology richer. 
If the cooling mechanism is efficient, the gravitational energy of the accreting gas is radiated away 
and the gas temperature is much smaller than the local virial temperature. This case is similar to 
the one involving objects with a solid surface: rj is equal to the binding energy per unit mass at the 
Innermost Stable Circular Orbit (ISCO), since we presume that the gas inside the ISCO falls quickly 
into the BH and is unable to emit further radiation. So, for Schwarzschild BHs r/ = 0.057, while 
for Kerr BHs the efficiency parameter can be as high as 0.42 (Shapiro & Teukolsky (1983)). On the 
other hand, if the cooling is not efficient, the gravitational energy is stored as thermal energy into 
the gas rather than being radiated. That can occur if the gas density is very low and particles do 
not scatter each other very much, or in the opposite case, when the medium is optically thick and 
radiation is trapped, as happens for high accretion rate. Here, unlike neutron stars, BHs have an 
event horizon and the energy can be lost into the BH. rj turns out to be very small and the accretion 
luminosity can be low. The accretion rate of matter can thus be high. 

We will argue that the BHs at the center of the first stars may accrete at the Bondi rate, with the 
star adjusting to keep the luminosity equal to the Eddington value, corresponding to a small value 
of f] in Eq. |3T|) . With Bondi accretion, the BHs can swallow the star in a short time, becoming 
10-1000 Mq BHs. There is considerable discussion of BHs accreting material inside stars in the 
literature. We present here some of the possibilities for the evolution of these objects. In all cases, 
the end result is a 10-1000 Mq BH. In the case of radiative stars, we may follow Begelman (1978), 
where the author discusses the steady flow accretion onto a Schwarzschild BH of a non-relativistic gas 
where the radiation pressure at infinity is much larger than the particle pressure and the radiation- 
particle coupling is provided by the Thomson scattering. The medium is optically thick at all the 
relevant scales and radiation is transported by diffusion and convection. Here one finds a trapping 
surface at the radius 

R, = MmfIh^ (32) 
47rmp 

inside which the radiation is convected inward and swallowed by the BH faster than it can escape to 
infinity|!| If Rt is much larger than the Bondi radiu Bb = 2GM/v'^, then the radiation is effectively 
trapped, that is it is convected inwards faster than it can diffuse outwards. In our case, using eq (|10p . 

Rb UO'°g/ Vlg/cmV ^3-107 cm/sy ' ^ ' 

Given the typical densities and temperatures inside Pop. HI stars, this condition is verified, the 
process is essentially adiabatic and in principle the BH is capable of accreting at an arbitrary high 
rate. Since radiation is trapped, the luminosity produced by the accretion process can not exceed the 
Eddington value, and the radiative efficiency effectively adjusts in order to keep L LEdd (Begelman 
(1978)). As long as accretion is spherical, with zero angular momentum, the central PBH can accrete 
ad libitum, and eventually swallow the whole star. In the presence of limited angular momentum we 
can argue that as long as the accretion disk that forms is all contained within the trapping radius, 
then radiation remains trapped and the growth of the BH can continue (Volonteri & Rees (2005)). 
We can take as a safe limit the condition that the disc is all within the trapping radius; this provides 
a lower limit to when accretion stops. The outer edge of the accretion disk, is roughly where the 
specific angular momentum of the gas equals the angular momentum of a gas element in a Keplerian 



'Clearly Rt cannot be larger than the radius of the star, R*. In this case, we take Rt = iJ* 
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circular orbit, therefore: 



Rd 
Rb 




(34) 



where Cg is the sound speed and V{Rb) is the rotational component of the velocity at the Bondi 
radius. In this case it still seems possible that the radiative efficiency drops so that the BH can 
accept the material without greatly exceeding the Eddington luminosity. Relaxing the assumptions 
of zero angular momentum and absence of mechanical turbulence and/or magnetic fields, the actual 
matter accretion rate presumably decreases, but the evolution of the star is slowed down as well. 
On the other hand, for very high angular momenta, it sounds reasonable that the system looks like 
a collapsar (MacFadyen & Woosley (1999)). 

Convective stars: 100 Mq Pop. Ill stars are primarily convective (Heger et al. (2007)). One 
can compute the Eddington luminosity in the case of a BH inside a mostly convective star with a 
radiative outer envelope as follows (Begelman et al. (2008)). There is no radiation pressure inside 
the convective zone, so the luminosity from the BH can easily get to the radiative outer envelope. 
Out there radiation pressure does exist. Then the Eddington luminosity at this outer region (which 
basically contains the entire star) is the relevant quantity. In short, one should use the Eddington 
luminosity of the star rather than Eddington luminosity of the BH, which means substituting Af, 
for Mbh in Eq. ([8]). Doing this, one finds 



This value is significantly larger than the numbers obtained in Eq. ([8]) because it is the Eddington 
luminosity of the star rather than that of the BH. Here the accretion luminosity is bigger than the 
fusion luminosity. The consequence for the star will be that it must expand, it will cool, and fusion 
will shut off. At that point the star looks like the quasistars in Begelman et al. (2008). These 
authors have worked out the stellar structure for a BH of arbitrary mass inside a star of arbitrary 
mass, where the only heat source is accretion luminosity. These authors were studying a different 
problem: they were not looking at Pop. Ill stars in 10^ Mq haloes; instead they were looking at 
cooler regions of similar content in 10^ Mq haloes. Although the context was different, the resultant 
objects should be very similar. 

There are 2 possibilities for the accretion: 1) The accretion may be spherical. In that case 
rj can be very small, as discussed in Bisnovatyi-Kogan & Lovelace (2002). There is no problem 
having t] — 10~^ so that the Eddington luminosity in Eq. ([35]) is compatible with Bondi accretion 
at M — erg/sec. Then it takes a thousand years to swallow the 100 Mq star (see Eq. (53) 
in Begelman et al. (2008)). Even more interesting is to contemplate the possibility that the star 
is accreting further mass from the halo outside it, e.g. at a rate 0.01 Mq/ year (McKee & Tan 

(2007) Then the BH can end up very large as seen in Eq. (52) of Begelman et al. (2008), possibly 
eating all 10^ Mq of baryons in the DM halo. 

2) The accretion may be in a disk. If the disk is thin and radiatively efficient, then t] ~ 0.1 
and M <C Mb (the accretion rate is much slower than Bondi). However, in different geometries, ij 
can become much smaller (Abramowicz & Lasota (1980)). Begelman et al. (2008) claim that the 
accretion stops once you hit "the opacity crisis." This happens when the photosphcric temperature 
(at the edge of the star) goes down to a critical value, so that the radiation pressure in the outer 
envelope vanishes, nothing prevents the star from going super-Eddington and blowing off all its gas. 
This leaves behind an exposed BH that no longer accretes anything. They find that for a fixed stellar 
mass of 100 Mq, the resultant object is a 10 solar mass BH in 10^ years, but nothing bigger, due to 
this opacity crisis. On the other hand, if the star is accreting further material from the outside, then 
you can end up with a 400 Mq BH if the accretion rate of material onto the star is 10~^ Mq/ yr 
(McKee & Tan (2007)) or 4000 Mq BH if the accretion rate onto the star is 10"^ Mq/ yr. Again, it 
takes 10^ years to reach this. In the meantime, during this lO'^ years, you have a "PBH Dark (matter 

'^The accretion rate for Pop. Ill stars is still highly uncertain, and certainly variable as a function of time. The 
values we quote are higher than typical estimates for prolonged accretion rate (see e.g. Figure 8 of McKee & Tan 

(2008) ), but still definitely possible, especially if PBHs somewhat reduce feedback effects. 



Lbh = Le.* = 10*° erg/sec (A/,/100Mq) . 



(35) 
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powered) Star" , i.e. a Pop. Ill star powered by accretion luminosity rather than by fusion. The 
exact accretion rate is none the less quite uncertain. Convective energy transport is itself limited 
and bounds the accretion rate (Begelman et al. (2008)) 

Mbh < . (36) 

Since Cg ~ 10~^ — 10~^, the actual accretion rate might be much smaller than the Bondi rate Mb, 
unless r] is quite small, say rj < 10"'* — 10~^. This is not a problem for spherical accretion, but might 
affect results for disk accretion. Regardless, this will require more study. Even accretion onto the 
first stars without the additional effects from primordial black holes is presently still an unsolved 
problem. 

We have argued that the BHs at the center of the first stars may accrete at the Bondi rate, so 
that the BHs can swallow the star in a short time, becoming 10-1000 Mq BHs. This mechanism 
may produce the seeds to generate the super-massive BHs which have been observed even at high 
redshifts and at the centers of galaxies. 



6 Summary and conclusions 

Primordial black holes in the mass range Mpbh ~ 10*^ — 10^^ g are viable dark matter candidates. 
They may be produced in the early Universe by many mechanisms and so far there are no constraints 
on their possible abundance. Assuming that they make part of the cosmological dark matter, wc 
expect that due to dynamical friction primordial black holes will make up a small but significant mass 
fraction of the first stars. Primordial black holes with masses smaller than about 10^^ g do not have 
a significant effect on the evolution of primordial stars, because their timcscalcs for Bondi accretion 
and for dynamical friction are larger than the lifetime of a Main Sequence star of 10 — 100 Mq. 
On the contrary, primordial black holes heavier than 10^^ g might sink quickly to the center of the 
star by dynamical friction and form a larger black hole, which could swallow the whole star in a 
short time. So, Pop. HI stars would likely have lived for a short time, with implications for the 
reionization of the Universe after the cosmic dark ages and the nature of the first supernovae; in 
fact they may preclude any supernovae from the first stars. Although the BH swallowing the star 
shortens the star's lifetime and its contribution to reionization, the newly formed hole can become a 
new, alternative source of ionizing photons. The 10-100 Mq BHs that form by swallowing the Pop. 
Ill stars may grow even larger: they reside in 1000 Mq of gas that are in excess of the Jeans mass 
and may fall into the BH. Black holes of mass 1-1000 Mq may result. 

Depending on the accretion mechanism at this point, the black hole may accrete more matter 
and grow larger. The 10^ Mq minihalocs of dark matter contain ^ 10^ Mq of baryonic matter. This 
accretion from the minihalo, as well as from other haloes merging with the one containing the black 
hole, would be from low density gaseous material {p ~ 10~^^ g/cm^), which is considerably different 
from the accretion we considered earlier from within the star (p ^ I g/cm^). In the case of accretion 
from the low density gas outside the star, feedback may become important. As we have shown, the 
timescale for the Pop. Ill stars to become black hole can be much shorter than the lifetime of the 
Pop. Ill stars (3 Myr for a 100 Mq star), so that the feedback due to stellar heating and ionization 
of the medium surrounding the black hole may be minimal. However, the accretion may well be in 
a disk, with the accompanying radiation pressure as well as radiative feedback due to the accretion 
(Silk & Rces (1998); Springcl et al. (2005); Giotto & Ostrikcr (2001); Li et al. (2007); Pclupcssy 
et al. (2007); Alvarez et al. (2008)) limiting the accretion rate. A recent study (Alvarez et al. 
(2008)) of the radiative feedback from the black hole accretion has been done for the case off] = 0.1 
and a Pop. Ill star that has imdergone its full lifespan, and finds reduced accretion onto the BH; 
the story may be different here. We have not studied these later stages. Since the end-products are 
10 — 10^ Mq black holes, these objects may serve as seeds for Intermediate Mass Black Holes; the 
super-massive black holes which have been seen already at high redshifts (Haiman & Loeb (2001); 
Volonteri & Rees (2006)) and may be the progenitors of the super-massive black holes which are in 
the center of every normal galaxy today. 
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Even if the primordial black holes do not explain the entirety of the dark matter in the Universe, 
they may still play a role in the first stars. Heavier primordial black holes than the ones studied 
here, i.e., primordial black holes with Mbh > 10^^ g, are observationally constrained to be only 
a fraction of the total dark matter in the Universe, and yet could be important in the first stars. 
It would only take one such black hole to be pulled into the star via dynamical friction (timescale 
~ 10^ yr for a 1 Mq black hole to get from 1 pc out into the center of the star (see Eq. ((30| ) 
and to quickly eat up the whole star. In fact, a single massive primordial black hole would already 
have a major effect during the protostellar phase while the molecular cloud is collapsing down into 
a protostar: the molecular cloud would already collapse into a black hole. In this case the fusion 
phase of a Pop. Ill star would be completely avoided. Another possibility would be to have the 
dark matter consist primarily of Weakly Interacting Massive Particles (WIMPs) but with a small 
component of primordial black holes. In that case there would be dark stars powered by WIMP 
annihilation (Spolyar et al. (2008)), which would become black holes once the primordial black holes 
sink to the center of the dark star. 

In principle, if the effects described in this paper do not take place, one could place bounds on the 
black hole abundances of various masses. For example, if primordial black holes swallowed primordial 
stars too quickly, the cosmological metal enrichment would be problematic and in absence of viable 
alternatives, the current allowed mass range Mpbh ^ 10^^ — 10^^ g could be further reduced to 

_ 1017 _ 1022 g. 
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Radius (cm) 

Figure 1: Density (top panel) and temperature (bottom panel) profile for the n ~ 3 polytrope star 
of mass 100 Mq used in our simulations. 
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